A mean field study of single-particle spectra evolution in Z = 14 and N = 28 chains 
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We study the mechanisms which reduce the proton 1^3/2-1^5/2 spin-orbit splitting and the neutron 
I/7/2 subshell closure in 42 Si. We use various self-consistent mean field models: non-relativistic 
Skyrme-Hartree-Fock and relativistic density-dependent Hartree-Fock. Special attention is devoted 
to the influence of a tensor component in the effective interaction. It is found that the tensor force 
indeed governs the reduction of the Id proton spin-orbit splitting. On the other hand, the reduction 
of the neutron I/7/2 subshell closure is not clearly related to the tensor force. 

PACS numbers: ... 



I. INTRODUCTION 



With the progress in experimental studies of single- 
particle nuclear spectra it has become clear that the 
single-particle levels undergo modifications when the 
neutron number TV or proton number Z change. From 
the experimental side, these studies are a difficult task 
because the nuclei involved are generally short-lived and 
their production rate can be low. It is important to 
have in parallel theoretical models which may serve as 
guidelines for exploring the evolution of nuclear proper- 
ties with N and Z numbers. A recent example is the work 
of Ref. [l[ drawing the attention on the possible effects of 
the neutron-proton tensor force on single-particle spectra 
along an isotopic or isotonic chain. This has triggered a 
number of experimental 0, 0] and theoretical [Jd, @, 0] 
works exploring this issue. 

In the present work we would like to use the self- 
consistent mean field approach to analyze two findings 
based on experimental data complemented by a shell 
model analysis @, Q: 1) the shrinkage of the proton 
spin-orbit splitting 1^3/2-1^5/2 in the Z = 14 isotopes 
when going from N = 20 ( 34 Si) to N = 28 ( 42 Si); 2) 
the quenching of the neutron shell closure gap 2p 3 / 2 - 
I/7/2 when going along the isotonic chain N = 28 from 
Z = 20 ( 48 Ca) to Z = 14 ( 42 Si). Our analysis is 
based on both non-relativistic Hartree-Fock-BCS (HF- 
BCS) approach with Skyrme-type interactions [|[ and 
relativistic Hartree-Fock-BCS (RHF-BCS) approach 0. 
Our aim is to discuss the effects, among other things, of 
the tensor component of the two-body effective interac- 
tion. The HF-BCS model is well suited for that purpose 
with the several Skyrme-type tensor forces proposed re- 
cently. On the relativistic side, the widely used relativis- 
tic mean field (RMF) model cannot help because it is 
just a Hartree approximation where a tensor component 
in the interaction would have no effect. Therefore, we use 
in our analysis the RHF model where exchange (Fock) 
terms are kept and the pion-induced tensor interaction 
can play its role [lj| • 



This paper is organized as follows. In Sec. II we briefly 
give the main ingredients of the two models used in this 
analysis, namely the non-relativistic HF-BCS model with 
a Skyrme type interaction, and the relativistic RHF-BCS 
model with density-dependent effective Lagrangians. In 
Sec. Ill the results obtained for the Z = 14 isotopic chain 
and the N — 28 isotonic chain of nuclei are discussed. 
Conclusions are drawn in Sec. IV. 



II. SELF-CONSISTENT MEAN FIELD MODELS 

We assume a spherical symmetry description of the 
nuclei studied here. Although this is a strong assumption 
for some nuclei of the isotopic and isotonic chains that 
we consider, it allows an easier insight into the role of 
separate parts of the mean field such as the central and 
spin-orbit potentials. 



A. HF-BCS mean field with Skyrme interaction 

The Skyrme-HF model is widely used and we refer to 
Ref. [1] for the notations of the present work as well as 
for the effective interaction SLy5 used here for generating 
the self-consistent mean field. In addition to the usual 
central and spin-orbit components of the Skyrme interac- 
tions (represented here by the SLy5 parametrization) we 
want to study the possible effects of a tensor component 
of the Skyrme force. This component was introduced in 
earlier versions of the Skyrme force [ll|, [l2|, [l3[ but it is 
only recently that attention was focused on its effects on 
single-particle spectra in spin-unsaturated nuclei. 

Along the chain of N = 28 isotones (Z = 14 isotopes) 
the proton-proton (neutron-neutron) pairing correlations 
will play a role, and they are treated in BCS approxima- 
tion [lJI using a zero-range, density-dependent pairing 
interaction of the form: 



y(n or p) _ y(n or p) / ^ 



p(R) 

Po 



<Kn-r 2 ) , (i) 
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where R 



(n + r 2 )/2, po 



0.16 fm 3 is the nu- 



clear matter saturation density and V^ n or ^ is fitted for 
each chain to reproduce the empirical two-neutron(two- 
proton) separation energies. The pairing window con- 
tains the orbitals of the lp, 2s-\d and 2p-lf shells. 

The effect of a tensor component in the Skyrme in- 
teraction is very simply described by a modification of 
the coefficients a and (3 which multiply neutron or pro- 
ton spin densities as explained, e.g., in Ref. This 
modification affects the single-particle spin-orbit poten- 
tials, especially in the case of spin-unsaturated subshclls 
which contribute importantly to the spin densities. One 
has: 



a = ac + &t , 
P = Pc + Pt , 



(2) 



where otc,Pc come from the central, velocity-dependent 
part of the Skyrme force whereas or, Pt come from the 
tensor component. Using the Skyrme force SLy5 the val- 
ues of a-r and were determined Q to be ax — —170 
McV fm 5 , P T = 100 MeV fm 5 , and we shall adopt these 
values. 

In the non relativistic Skyrme-HF model the radial 
HF equations can be expressed in terms of an energy- 
dependent equivalent potential V^: 



2m 



1(1 + 1) 



Mr) 



■V«(r,e)^(r)=e^(r) , 
(3) 



where 



V£(r, e) = ^U (r) + ^E#(r) + V™ . (4) 
M m m M 

Here, Ug(r) = U so (r) x + 1) - 1(1 + 1) - 3/4], U so (r) 
is the spin-orbit HF potential, 



B. RHF-BCS mean field with density-dependent 
Lagrangians 

Another type of mean field approach is based on covari- 
ant effective Lagrangians. A very popular and success- 
ful version is the RMF (see, e.g., Ref. [HI). However, 
the RMF is only a Hartree approximation and it can- 
not describe the effects of a two-body tensor interaction 
brought about by exchanging pions, or rho mesons with 
tensor coupling because these effects would appear in the 
exchange (Fock) contributions to the total energy. Thus, 
the RMF aproach is not suitable for our present purpose. 

Recently, progress have been made with the RHF-BCS 
approach. New density-dependent effective Lagrangians 
have been proposed and they are able to give a satis- 
factory overall description of bulk properties of nuclei 
throughout the mass table 0, [l7| • In this work we will 
compare results obtained with two different effective La- 
grangians: the parameter set PKOl 0| which includes 
a, to, p- vector and tt exchanges, and the parameter set 
PK02 which does not include 7r exchange. Thus, the 
comparison of the predictions of PKOl and PK02 can 
give an idea of the effects on single-particle spectra due 
to a tensor component in the interaction, since the Fock 
terms of the one-pion exchange are indeed sensitive to 
this tensor component. 

The two effective Lagrangians PKOl and PK02 are 
treated in the RHF-BCS approximation within a spheri- 
cal symmetry assumption as in the Skyrme HF-BCS case. 
The pairing interaction has the same functional form as 
in Eq.JT]), and the interaction parameters V^ n ° r p ^ are 
adjusted as in the non-rclativistic case. The RHF-BCS 
equations are self-consistently solved in coordinate space 

u. 



III. RESULTS AND DISCUSSION 
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and Uo(r) and m*(r) are the HF central potential and 
effective mass, respectively The Coulomb potential 
is included in Uo, in the case of protons. With the help of 
Eqs. ([3j)- (0) we can write e as a sum of kinetic, central, 
spin-orbit and non-locality (or m*-dependent) contribu- 
tions: 



Cki: 



(6) 



In the above expression the three last terms correspond 
to the three terms of Eq. ([4]) . This decomposition is useful 
for understanding the evolution of single-particle energies 
when N or Z is changing [l5[ , and it will be used in the 
discussion of results in Sec. III. 



A. The 2s-ld proton spectra in Z = 14 isotopes 

We first examine the evolution of the proton single- 
particle spectra of the 2s-ld shell in the isotopic chain 
Z = 14 when the neutron number increases from N = 
20 ( 34 Si, ulf r /2 empty) to TV = 28 ( 42 Si, ulf 7/2 filled). 
The question of reduction of the proton 1^3/2-1^5/2 spin- 
orbit splitting when going from 34 Si to 42 Si is particularly 
interesting. In a shell model analysis of the data it is 
found @ that this spin-orbit splitting is reduced by about 
1.94 MeV if one fills the I/7/2 neutron subshell. 

In the left panel of Fig. [T] we present the 2s-1g? proton 
levels calculated with the Skyrme force SLy5, with and 
without the tensor component of Ref. [i| . Corresponding 
results obtained with the RHF-BCS model are shown in 
the right panel, using the PKOl (with pion exchange, i.e., 
with tensor coupling) and PK02 (without pion exchange) 
paramctrizations. The general down-going trend of all 
results as a function of N is easily understood as an effect 
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FIG. 1: The 2s-ld proton single-particle levels in the Si isotopic chain calculated in the mean field approach. Left panel: 
non-relativistic model with SLy5; right panel: relativistic model with PKOl and PK02. 



of the neutron-proton symmetry potential. 



TABLE I: The ld 3/2 -ld 5 / 2 splitting for for the nuclei a4 - 42 Si 
calculated for the models described in the text; The units are 
MeV. 



(N,Z) 


(20,14) 


(22,14) 


(24,14) 


(26,14) 


(28,14) 


SLy5 


6.38 


6.33 


6.35 


6.35 


6.36 


SLy5+T 


8.62 


8.26 


7.85 


7.43 


7.10 


PKOl 


7.26 


6.81 


6.41 


6.07 


5.83 


PK02 


7.46 


7.24 


7.02 


6.82 


6.65 



The calculated values of the proton spin-orbit splitting, 
A so (7V) = e ld 

3/2 ~ e id 5/2 \N are shown in Table [J It can 
be seen that the reduction of A so (7V) from N = 20 to 28 
is practically zero for SLy5 but it becomes 1.52 MeV for 
SLy5 plus tensor. For PKOl and PK02 the reduction 
is 1.43 MeV and 0.81 MeV, respectively. It seems that 
the tensor force helps to reduce A so (N) in both the non- 
relativistic and relativistic models. 

Actually, the two relativistic models PKOl and PK02 
differ not only by the pion-coupling term but also by the 
strengths of other meson-nucleon couplings and there- 
fore, one cannot relate the changes in A so (N) only to 
the tensor coupling. On the other hand, it is an eas- 
ier task to analyze the evolution of A so (N) in the non- 
relativistic model. We use the decomposition of Eq. ^ to 
cast A so (N) into a sum of 3 terms: 1) a kinetic plus cen- 
tral contribution; 2) a spin-orbit contribution; 3) a term 
due to m*(r) < m and related to the non-locality of the 
mean field. These three terms are shown in Fig. O The 



kinetic and central contributions are of opposite signs and 
comparable magnitudes, their sum remains small from 
A = 20 to 28 and the tensor force has little effect on 
it. The spin-orbit term is nearly flat for SLy5 whereas 
the results of SLy5 + tensor show a strong decrease of 
about 1 MeV. This effect is easily understood because of 
the properties of the Skyrme-HF mean field and of the 
neutron-proton tensor interaction. Firstly, the neutron 
spin density J„ becomes larger when one fills the neutron 
orbital I/7/2 [8|- Then, the proton spin-orbit potential 
which contains a term PtJu changes accordingly. Since 
the neutron I/7/2 orbital is of the type j = I + 1/2, the 
sign of this change is such that the proton ld 5 / 2 state 
is pushed up and 1^3/2 is pushed down according to the 
general rule of matrix elements of a tensor interaction [l|] . 
The proton spin density remains the same for the entire 
chain, so the term axJ P contributes in a same way to all 
nuclei. Finally, the non-locality term is flat for the SLy5 
results but it decreases by about 300 keV from A^ = 20 
to 28 for the SLy5 + tensor case because it is practically 
proportional to A so (N). 

Thus it can be concluded that, in the framework of the 
self-consistent mean field with Skyrme interactions, the 
reduction of the proton 1^3/2-1^5/2 spin-orbit splitting 
in 42 Si can be attributed to the effect of a tensor com- 
ponent in the effective interaction. For the RHF-BCS 
model, one observes that the reduction is also enhanced 
by a tensor component brought about by a pion-induced 
interaction. However, all of the present discussion is done 
strictly in the framework of a static mean field approach, 
and it must be kept in mind that effects beyond mean 
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We observe again the down-going trends as functions of 
Z, a manifestation of the increasing symmetry poten- 
tial in the mean fields. Looking at the N = 28 gap 

A ga p(Z) = £2 P3/2 - ei/r/aU we can see n ' om 3 that, 
in the non-relativistic case adding a tensor force affects 
little the 2p 3 / 2 state but pulls down the I/7/2 state and 
consequently increases A gap (Z) for all Z values. On the 
other hand, in the relativistic case the position of the 
I/7/2 state is almost the same for both PKOl and PK02 
parametrizations but the 2p 3 / 2 state is lower with PKOl 
(with tensor component) than with PK02, so that the 
model with a tensor component seems to lead to smaller 
Agap(^). However, as noticed above, the differences be- 
tween PKOl and PK02 arc not only in the tensor com- 
ponent contrarily to SLy5 and SLy5 + tensor. 



TABLE II: The 2p s/2 -lf 7 /2 gap for the isotonic chain N = 28; 
The units are MeV. 



(N,Z) 



(28,14) 



(28,16) 



(28,18) 



(28,20) 



SLy5 
SLy5+T 
PKOl 
PK02 



2.55 
3.79 
2.13 
2.93 



3.09 
4.08 
3.52 
4.16 



3.49 
4.66 
5.11 
5.26 



3.75 
5.15 
4.78 
5.22 



FIG. 2: The contributions to A so (iV) according to the de- 
composition Eq.©. The results correspond to the Skyrme 
interaction SLy5, with (solid lines) and without (dashed lines) 
tensor component. 



field like particle- vibration coupling [18| have to be eval- 
uated for these nuclei. Because of these missing effects 
one cannot compare directly our calculated spectra with 
empirical ones. Nevertheless, the predictions of the mean 
field approach concerning evolutions of energy differences 
remain meaningful. 



B. The 2p-l/ neutron spectra in N = 28 isotones 

We now discuss the changes of neutron single-particle 
energies in the 2p-\f shell when going along the TV = 28 
isotones from 48 Ca (Z = 20) to 42 Si (Z = 14). 

Experimental observations combined with a shell 
model analysis @, Q lead to the conclusion that a com- 
pression of the neutron 2p-lf spectrum occurs when one 
removes protons from the W3/2 and 2si/ 2 orbitals, keep- 
ing N = 28 fixed. This compression could affect by about 
1 MeV the distance between the neutron I/7/2 and 2p 3 / 2 
states and reduce the N = 28 gap from 4.8 MeV in 48 Ca 
to about 3.8 MeV in 42 Si. 

In Fig. [3] we present the energies of the three low- 
est states of the neutron 2p-\f shell calculated with 
the non-relativistic and relativistic mean field models. 



The calculated values of the N — 28 gap A gap (Z) are 
displayed in Table [ill It can be seen that all models pre- 
dict a substantial reduction of this gap when going from 
Z = 20 ( 48 Ca) to Z = 14 ( 42 Si). The largest reduc- 
tion is for the two models PKOl (2.65 MeV) and PK02 
(2.29 MeV), whereas the Skyrme models SLy5 and SLy5 
+ tensor predict a reduction of 1.20 MeV and 1.36 MeV, 
respectively. These values are consistent with the exper- 
imental estimate of 1 MeV. The non-relativistic results 
also indicate that the tensor component of the interac- 
tion does not play an essential role in this reduction. This 
may be supported by the fact that both PKOl and PK02 
relativistic models lead to very similar reduction. 

One can have a closer look at the mechanism of this 
gap reduction by performing for the non-relativistic mod- 
els an analysis similar to that of Subsec.III.A. In Fig. 0] 
are shown the contributions to A gap (Z) coming from the 
different terms of Eq. © . Comparing the results of SLy5 
and SLy5 + tensor one can see that the tensor force has 
an overall effect on the spin-orbit contribution to A gap (Z) 
but it affects little the gap reduction from Z = 20 to 
Z = 14. The other two terms practically determine the 
gap reduction. In the SLy5 + tensor case for example, 
the partial decrease from Z = 20 to Z = 14 due to the 
kinetic + central effect is 840 keV, while the non-locality 
effect is 564 keV. 

We stress again that, in the present study the particle- 
vibration effects are still missing and their effects on 
Ag a p(Z) remain to be evaluated. 
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FIG. 3: The 2p 1( / 2 , 2p 3 / 2 and I/7/2 neutron single-particle levels in the N = 28 isotonic chain calculated in the mean field 
approach. Left panel: non-relativistic model with SLy5; right panel: relativistic model with PKOl and PK02. 



IV. CONCLUSION 

In this work we have studied the mechanisms which 
lead to a reduction of the proton 1^3/2-1^5/2 spin-orbit 
splitting A so (JV) with decreasing the mass number in the 
Si isotopic chain. We also investigate a quenching of 
the neutron 2p 3 / 2 -l/7/2 g a P A gap (Z) as one goes from 
48 Ca towards 42 Si. We have used the self-consistent 
mean field approach in its non-relativistic version with 
the SLy5 parametrization, and its relativistic covariant 
version with exchange (Fock) terms. 

One of the goals of this study is to determine to 
which extent a tensor component in the effective nucleon- 
nucleon interaction can affect the values of A so (iV) and 
Ag a p(if). Our main conclusions arc: 1) the reduction of 
A so (iV) when going from 34 Si to 42 Si is mainly due to 
the presence of a tensor component in the effective inter- 
action. The magnitude of the change in A so (A) is con- 
sistent with the empirical observations. 2) On the other 
hand, the evolution of A gap (Z) when Z decreases from 20 
to 14 does not depend strongly on a tensor component in 
the interaction. These conclusions can be reached by us- 
ing either a non-relativistic Skyrme mean field approach, 
or a relativistic Hartrcc-Fock framework. 



In the case of the Skyrme- Hartree- Fock model, we have 
performed a detailed analysis of the origin of the evolu- 
tion of A gap (Z) and A so (N). The reduction of the former 
can be traced back to changes occurring in the symme- 
try part of the central potential, to the m*-term of the 
Skyrmc-HF model, and only to a small extent to the 
spin-orbit potential. 

It would be important to extend in future studies 
this type of analysis to include effects such as particle- 
vibration coupling which are beyond the present mean 
field models and which are known to affect single-particle 
spectra,. 
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FIG. 4: The 3 contributions to A gap (Z) according to the 
decomposition Eq.©. The results correspond to the Skyrme 
interaction SLy5, with (solid lines) and without (dashed lines) 
tensor component. 



